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Abstract. The imprint of neutrino oscillations on the Cosmic Microwave Background (CMB) anisotropy and po- 
C^l ' larization power spectra is evaluated in a ACHDM model with two active neutrino flavors, consistent with the 

structure formation models and the atmospheric neutrino oscillations data. 

By using the Fisher information matrix method we find that the neutrino oscillations could be detected by Planck 
tH- , surveyor if the present value of the lepton asymmetry L^ > 4.5 x lO^'^, the difference of the neutrino squared masses 

OS ' Ato^ > 9.7 X 10~''eV^ and the vacuum mixing angle sin^ 26*0 > 0.13, showing the existence of a significant overlap 

^D , between the region of the oscillation parameter space that can be measured by Planck surveyor and that implied 

^^ ' by the atmospheric neutrino oscillations data. 

^. . . 
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1. Introduction 

•i-H . 

j^ ' The atmospheric neutrino results from Super-Kamiokande (Fukuda et al. 1998) and MACRO (Ambrosio et al. 

JH I 1998) experiments indicate that neutrinos oscillate. Those data are consistent with t/^ <-> i>r oscillations, but do 

not exclude i/^ oscillations to a sterile neutrino I's (see e.g. Foot et al. 1996; Foot & Volkas 1997,1999). The small 
value of the difference of the squared masses (5 x 10~'*eV < Am^ < 6 x lQ~^eV ) and the strong mixing angle 
(sin^ 29 > 0.82) suggest that these neutrinos are nearly equal in mass as predicted by many models of particle 
physics beyond the standard model (see e.g. Primack & Gross 1998 and the references therein). Also, the LSND 
experiment (Athanassopoulos et al. 1998) support I'f^ ^-> Ve oscillations with Am^ < 0.2eV^ and other different 
types of solar neutrino experiments (Bahcall et al. 1998) suggest that z/g could oscillate to a sterile neutrino I'e <-^ i^s 
with Am^ ~ lO-^eV^. 

The direct implication of neutrino oscillations is the existence of non-zero neutrino masses in the eV range, and 
consequently a not negligible hot dark matter contribution to the total mass density of the universe (i.e. a density 
parameter fl^ 7^ 0). 

In view of the uncertainty in the direct neutrino experiments, a possible significant contribution can be obtained 
from the study of the cosmological implications of neutrino oscillations. In particular, the Cosmic Microwave 
Background (CMB) anisotropy pattern reflects the conditions in the universe at the time of the last scattering 
between photons and electrons, occurring at a redshift z « 1000 for standard recombination models. This means 
that all the physical processes occurring before this epoch could have left imprints on the CMB angular power 
spectra. 



** Present Address: Physics Dept., Purdue University, West Lafayette, IN 47907-1396, USA 
Send offprint requests to: lpopa@venus.nipne.ro or popa@tesre.bo.cnr.it 



2 L. Popa et al.: On the detection of neutrino oscillations with Planck surveyor 

The standard Cold Dark Matter (CDM) model normalized to COBE/DMR data (Smoot et al. 1992, Wright 
et al. 1994, Gorski et al. 1994, Bennet et al. 1996a) predicts both the amplitude and the shape of the CMB power 
spectrum at small scales inconsistent with the observations of the Large Scale Structure (LLS) of the universe 
as derived by galaxy surveys (e.g. Scott & White 1994; White et al. 1995; Primack et al. 1995). Recent works 
(Primack 1998; Gawiser & Silk 1998) show that the Cold + Hot Dark Matter (CHDM) model is the single model 
with rim — 1 whose predictions agree with both LSS observations and current CMB anisotropics data. The C^^^DM 
model (Primack et al. 1995) with two 2.4eV neutrinos contributing with 20% to the total density of the universe 
and with another 80% contribution from cold dark matter and a small baryon fraction, agrees remarkably with all 
available observations only if the Hubble parameter is Hq = 50 km s^^Mpc^^ {ho ~ iJo/100 km s~^Mpc^^ = 0.5). 

However, increasing evidence for larger value of the Hubble parameter, ha ~ 0.6 — 0.8, (see e.g. Fukugita, Liu 
& Sugiyama 1999 and the references therein) are inconsistent with standard CHDM models. Such high values of 
the Hubble parameter make an Q„i = 1 universe suspiciously young unless the cosmological constant is non-zero. 
A way to improve the agreement of the standard CHDM models with eV neutrino mass with high values of Hubble 
parameter is to consider the relic neutrino degeneracy (the degeneracy parameter is defined as: ^^ — ^^/T^,^ where 
/Xjy is the neutrino chemical potential and T^ is the neutrino temperature), that enhances the contribution of 
neutrinos to the total energy density of the universe (Larsen & Madsen 1995). The cosmological implications of the 
neutrino degeneracy has been often considered in the literature: it can change the neutrino decoupling temperature 
(Freese et al. 1983; Kang & Steigman 1992), the abundances of light elements at the big bang nucleosynthesis (BBN) 
(Steigman et al. 1977; Kang & Steigman 1992), the CMB anisotropics and the matter power spectrum (Kinney & 
Riotto 1999). For the CHDM models with degenerated neutrinos, the increase of the neutrino chemical potential 
increases the effective number of neutrino species in the relativistic era (Larsen & Madsen 1995, Hannestad 2000), 
and can bring the power spectral shape parameter in the range required by observations [for CDM models with 
low baryonic content the power spectral shape parameter (which basically measures the horizon scale at matter- 
radiation equality) is defined as (Dodelson et al. 1994): F « i7m/io((7*/3.36)~^/^, where g* counts the relativistic 
degrees of freedom and g, = 3.36 corresponds to the standard model with photons and three massless neutrino 
species; observations require 0.22 < F < 0.29]. The free-streaming properties of degenerated neutrinos also differ: 
when the neutrino chemical potential is increased, the neutrino number density is higher, although the neutrino 
energy density and the mean momentum change a little. The effect is the increase of the radiation energy density 
and the delay of the matter-radiation equality. It is shown (Larsen & Madsen 1995) that an Q.m = 1 CHDM model 
with one 2.4eV neutrino and ho — 0.7 is in good agreement with the observed LSS power spectrum of density 
fluctuations if the neutrino degeneracy parameter is ^^ w 2.8. Constraints on neutrino degeneracy coming from 
BBN (Kang & Steigman 1992) indicate -0.06 < ^^^ < 1.1 and \£,y^^ < 6.9. Recent works (Pal & Kar 1999; 
Lesgourgues & Pastor 1999; Kinney & Riotto 1999, Lesgourgues et al. 2000) derive bounds on neutrino degeneracy 
parameter in agreement with BBN predictions, by combining the current observations of the CMB anisotropics and 
LSS data. The likelihood analysis of the CMB anisotropy data obtained by the Boomerang experiment indicates 
bounds on massless neutrino degeneracy parameter (Hannestad 2000) of \^^^ ^ | < 3.7 if only one massless neutrino 
species is degenerated and \£,i,\ < 2.4 if the asymmetry is equally shared among three massless species. 
On the other hand, evidences has been accumulated that we live in a low matter density universe (see e.g. Fukugita, 
Liu & Sugiyama 1999 and the references therein). Indications like Hubble diagram of Type la supernovae (Riess et 
al. 1998, Perlmutter et al. 1998) and the acoustic peak distribution in the CMB anisotropy power spectra (Hancock 
et al. 1998, Efstathiou et al. 1999) point to a universe dominated by vacuum energy (cosmological constant A) that 
keeps the universe close to flat. The combined analysis of the latest CMB anisotropy data and Type la supernovae 
data (Efstathiou et al. 1999) indicates Qm — 0-25^012 ^^'^ ^A — 0-63iQ'23 (95% confidence errors) for matter 
and vacuum energy densities respectively. These values are close to those favoured by other arguments (see e.g. 
Efstathiou et al. 1999 and the references therein) like the ages of globular clusters, observations of large scale 
structure, baryon abundance in clusters. 

Adding a Hot Dark Matter component to the ACDM model (ACHDM) leads to a worse fit to LSS and CMB data, 
resulting in a limit on the total neutrino mass (Gawiser 2000) of rrii, < 2eV for a primordial scale invariant power 
spectrum and rrii, < 4eV for a primordial scale free power spectrum. A stronger upper limit is obtained (Fukugita, 
Liu & Sugiyama 1999) from the matching condition of the LSS power spectrum normalization as (defined as 
the rms amplitude of the galaxy power spectrum in a sphere of radius 8h^ Mpc) at the COBE scale and at the 
cluster scale. For the case of a ACHDM model having ilm — 0.3, Ha — 0.7, and a primordial scale invariant power 
spectrum, it is found an upper limit of the total non-degenerated neutrino mass ofrrii, < 0.6 if the Hubble constant 
Ho < 80 Km s-i Mpc^^ 
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In this paper we study the signature of rehc degenerated neutrino oscillations on the CMB anisotropy and 
polarization power spectra and address its detectability with the future CMB anisotropy space missions, MAP 
(Microwave Anisotropy Probe) (see Bennet et al. 1996b) and Planck (Mandolesi et al. 1998, Puget et al. 1998). 
We consider the impact of neutrino oscillations in the epoch after nucleosynthesis. Neutrino oscillations are mediated 
by weak interactions: this fact implies that the order of magnitude of the interaction rate is less than the Hubble 
expansion rate after nucleosynthesis. In order to get an appreciable effect we consider a large neutrino asymmetry 
as left by processes occurred before and during nucleosynthesis. An example of such large asymmetries is given 
by considering a relic neutrino degeneracy (Larsen and Madsen 1995), which has been recently subject of renewed 
interest (Kinney and Riotto 1999; Lesgourges and Pastor 1999; Pal and Kar 1999). We analyze the oscillations 
of relic degenerated neutrinos under the simple assumption of oscillations occurring between two active neutrino 
flavors 

Z^p ^ Vr i>^^ Vr 

which have the following mass hierarchy: 

We assume the third neutrino Vp as massless and non-degenerate. We consider a total neutrino mass contribution 
of m^ = vfii, + rtii,^ — 0.6 eV and a total neutrino degeneracy parameter ^^ = \^i, + i^^^\ < 5 accordingly to 
the present observational data summarized above. We assume a scale invariant primordial power spectrum, the 
presence of the scalar modes with spectral index n^ = 1 and ignore the contribution of the tensorial modes and 
the reionisation effects. This model is consistent with the large scale structure and CMB anisotropy data, allowing 
in the same time a pattern of neutrino masses consistent with the results from atmospheric neutrino oscillations 
experiments. 

In Section 2 we draw the basic formalism of the neutrino oscillation model necessary for understanding its impact on 
CMB angular power spectra. The features induced by the degenerated neutrino oscillations on the CMB anisotropy 
and polarization power spectra are presented in Section 3. In Section 4 we discuss the MAP and Planck capability 
of detecting neutrino oscillations. Finally, we summarize our main conclusions in Section 5. 
Throughout the paper we employ the system of units in which h = c = kB — i- 

2. Oscillations of degenerated neutrinos 

The neutrino oscillations take place due to the fact that neutrino mass eigenstate components propagate differently 
because they have different energies, momenta and masses. In the standard assumption of the mixing of massive 
neutrinos, the neutrino flavor eigenstates are described by a superposition of the mass eigenstate components 
(Particle Data Group 1998). For oscillations occurring in vacuum between two neutrino flavors v^ and Vt- the 
mixing can be written as: 

\i^^) = cos6'o|j^2) +sin6lo|i^3) (1) 

Wt) ^ -sin0o|'^2) +cos6'o|!^3), 

where 1^2 and V3 are the mass eigenstate components and 9a is the vacuum mixing angle. The mixing of antineutrinos 
can be obtained from equation (1) by performing the transformations v^ ^-s- Pfj, and Vr *-^ Vr- It is usual to consider 
that each neutrino/ antineutrino of a definite flavor is dominantly one mass eigenstate (Particle Data Group 1998). 
In this circumstance we refer to the dominant mass eigenstate component of Vy^jvy, as vi-, that of VtJVt as v-^ and 
to their difference of the squared masses as Am^ = m§ — vri^. 

The mixing of the mass eigenstate components are modified in the presence of the asymmetric background. The 
background mixing angle and the vacuum mixing angle are related through (see e.g. Foot et al. 1996): 

sin2 2e^ = — — , " ^ , . 2 

1 — 2zcos26'o + 2^ 

Here z = 2 <'pxVi,^ > /Arn^ (i = /x, r), < p > is the averaged neutrino momentum and < V^. > is the neutrino 
effective potential due to the interaction with the asymmetric background (Enqvist et al. 1992; Foot et al. 1996): 

< K. >= V2GFn^ (v - A,^^[n,^ + n^j") , (3) 
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where A^. is a numerical factor (A^^^- — 12.61), T^ is the neutrino temperature, Gp is the Fermi constant (Gp = 
1.17 X 10""Mey"2) Mz is the mass of Z boson, {Mz = 91.187GeV), n^ is the photon number density {n^ ~ 
T^/AA), n^. and riij^ are the neutrino and antineutrino number densities. 
The functions U arc defined as (Foot ct al. 1996; Foot & Volkas 1997): 

L' = i^, + L,,^ + L^^ + L,,^ +Le + v ii = Ai, r). (4) 

Here rj is approximately equal to the baryon to photon ratio, and L s are the lepton asymmetries defined as: 

La = , (5) 

where n s are the number densities. 

For the purpose of this work we consider large neutrino asymmetries (L^. > 10^''). Also, we consider that ly^ is 

non-degenerated and consequently L^^ = 0. In this case the functions L* may be written as: 

V~L,^+U^+L,^. (6) 

The effective potentials of antineutrinos can be obtained by replacing L* with — L* in equation (3). 

The neutrino flavor eigenstates oscillate via weak interaction processes. The average oscillation probability 
< V >ui of a flavor eigenstate |z-'i) {i — fj,,T) is defined as (see e.g.: Raffelt 1996, Foot et al.l996, Enqvist et al. 
1992): 

<V >ui^ sin^ 29,n / sin^ — ^ \ Tcouij^i -^ Vj) {i ^ fi,T j ^ ^,t i ^ j), (7) 

\ ^osc I 

where Vcouivi -^ Vj) is the neutrino averaged elastic collision rate, Lint is the mean distance between interactions 
(Foot ct al. 1996), and the neutrino oscillation length Lose is given by (see e.g. Raffelt 1996): 

Lose ^ -r-^ - 248cm 8 

A?7i^ MeV Am^ 

where E^. is the neutrino energy. We will consider the Lint » Lose and then (^sin^ {Lint /Lose)) -^ 1/2 (Foot et 
al. 1996). The averaged collision rate is given by (Enqvist et al. 1992): 

r(i., ^ i.,) = 0.77G|.T5l^, (9) 

no 

where hq is the number density of massless non-degenerated neutrino species. 

Equation (7) together with equations (2), (8), (9) gives the averaged oscillation probability < P >^. used in 
the next section to evolve in time the neutrino distribution functions. The oscillation probabilities of antineutrinos 
can be obtained by performing the transformations Vi ^^ Vi and Vj <-^ Dj in equations (7) — (9). 

3. CMB angular power spectra in the presence of neutrino oscillations 

The phase space distribution functions of neutrinos and antineutrinos /^. (i — v^^D^,VT^Dr) arc evolved indepen- 
dently, taking into account the redistribution of momenta due to oscillations and elastic collisions. 
We describe the time evolution of the phase space distribution functions in the expanding universe by the set of 
equations: 

^^^ = -<V>,^ U, {q)+ < V >.^ U^ (g), {i = /i, r; j = a^, r; i^ j) 



dt 



<v>-,^hM)+<v>u,hM)^ (10) 



where g = ap is the comoving momentum (Ma & Bertschinger 1995), p is the magnitude of the momentum 3-vector, 
a is the scale factor (a = 1 today), dt — da/{aH) and H is the expansion rate: 



^ ^ ^.Gptotja) ^ (11) 
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where: G is the gravitational constant and ptotio) is the total energy density. 

The equations (10) assume that no other physical processes take place apart from oscillations and elastic collisions, 

leading to the conservation of the total neutrino/antineutrino number densities {n^/ui^) and degeneracy parameters 

n-i, = Uy. + n^. , He = rip^ + n^^ , {i = /i, j = r) 

c. = e.. +c.,, 6 = 6, +&,, (12) 

where ^a ~ I^ua/Tu is the degeneracy parameter of the flavor a (a = fJ-iT) with the chemical potential /Iq, T^ is 
the neutrino/antineutrino temperature and |^p^| — |^!/^|. 

We have modified the CMBFAST code (Seljak & Zaldarriaga 1996) to take into account the evolution of the 
phase space distribution functions of neutrinos and antincutrinos at each value of the scale factor, by solving the 
set of equations (10). We start at the scale factor a ~ 10"^, when neutrinos and antincutrinos with the mass in eV 
range behave like relativistic particles with Fermi-Dirac phase space distributions: 

.U,i<l) = E IT -f — TT' •^'^'(9) = E IT +( — TT' (* " ^'^) (1^) 

where £,i is the degeneracy parameter, E^. = ^^q^ + a?ml. the neutrino/antineutrino energy, m^^ the dominant 
mass component of the flavor i, and q the comoving momentum. The system of neutrinos and antineutrinos drops 
out from thermal equilibrium with e*, photons, the small fraction of baryons and other massive species (e.g. /i^/i^ 
pairs) when the ratio of the averaged weak interaction rate to the expansion rate falls below unity. The decoupling 
temperature of degenerated neutrinos increases exponentially with the degeneracy parameter. For v^^ and v-r the 
decoupling temperature is given by (Freese et al. 1983, Kang & Steigmann 1992) Td ~ 10 MeV ^~2/3g(c/3)_ 
The present values of neutrino/antineutrino temperature T^,^^ and photons temperature Tq = 2.725±0.002K (Mather 
et al. 1999) are related through T^g/T^ ~ {S.9/g^,D)^ where g^,£! is the number of degrees of freedom in equilibrium 
at neutrino/antineutrino decoupling. For a degeneracy parameter ^^ < 15, Tjj ~ few MeV and the total number 
of relativistic degrees of freedom is g, ~ 43/4 (Freese et al. 1983). After the neutrino decoupling, the e+e^ 
pairs annihilate, heating the photons but not the neutrinos. From the entropy conservation one obtains that 
Ti/g ~ (4/ll)^/^To and the total degeneracy parameter £,u = £,v + iv^ is considered constant. 
We compute the mean energy density and pressure of each neutrino flavor as: 

Pvi+Pvi = TTT / dqq^E„.{fuM)+ hAl)), (* = A*,t) 



27r2 





4 fOC 



the neutrino and antincutrino number densities 

n4 fCO rp4: 



(14) 



n. 



■i ^ ^ d(lq^U^{q)l '^^' ^ 2n^ dqq'^fp^iq), (i = fj.,T) (15) 



and the present lepton asymmetry in the form of neutrinos (Kang & Steigmann 1992): 

i—f_i.T ^ \ / 



We also compute at the beginning of the code the present value of neutrino parameter density: 

t—fl.T 



n.= J2 ^^"^^"^°^\ (17) 



where Pc — 1.054 x 10^ /iq eV cm '^ is the critical density, and define the cold dark matter density parameter 

Following the same procedure used by the CMBFAST code, we compute in synchronous gauge the perturbations 
of the energy density, pressure, energy flux and shear stress (equations (52) — (56) from Ma & Bertschinger 1995) 
for neutrinos and antineutrinos of each flavor. 
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Fig. 1. The dependence of CMB anisotropy (left panels) and polarization (right panels) power spectra on Am^ ( panels al) 
and a2) ) sin^ 260 ( panels bl) and b2) ) and A^i, ( panels cl) and c2) ) for a ACHDM model having: ^6 = 0.023, ^a = 0.7, 
Q.C = 1 — 0,b — ^v — Ha, ho = 0.65, rriv ~ rrii,^ + m^^ = 0.6eV, ^^ = S,v^ + iv^ = 5. As a consistency test we present in 
panel al) the CMB anisotropy power spectra obtained for degenerate neutrinos with: rriv — 0.07eV, £,,^ — 3 (dashed line), 
nil, — 0.07eV, (,1, — (dash-dotted line) and ho = 0.65, fit = 0.05, Oa = 0.7, Q.c = \ — ^b — ^v — ^a. In panels bl) and 
b2) we present (dashed lines) the CMB angular power spectra obtained for: Q.^, ~ 0.015, ^v = 0, fib = 0.023, ^k = 0.7, 
fie = 0.262, ho = 0.65, one massless and two massive neutrino species (see also the text). We choose such values of neutrino 
oscillation parameters in order to clearly see the neutrino oscillation effects on CMB power spectrum. 



We performed the computations of the CMB anisotropy and polarization power spectra starting with Nq^^^ = 15 
equispaced sampling points of the phase space distribution functions and truncating the Boltzmann hierarchies 
for massive neutrinos (see equation (54) from Ma & Bertschinger 1995) at Imax — 15 for every q value. We then 
gradually increase these numbers until Nq^^^ = 100 and Imax = 50 for testing the numerical precision and obtaining 
a relative accuracy better than 10^'^ for the models with oscillations and without oscillations. 

Figure 1 presents the dependence of the CMB anisotropy and polarization power spectra on the neutrino 
difference of the squared masses Am^, the vacuum mixing angle sin^ 29o and the lepton asymmetry Li^, obtained in 
a ACHDM model having: fli, = 0.023, JIa = 0.7, J7c = 1 — fif, — Jl^ — f^A, /ig = 0.65. For all cases wc assume a total 
neutrino mass rrii, = m^^ + m^^ = 0.6eV and a total neutrino degeneracy parameter ^^ = ^^^ + ^^^ = 5. We also 
assume a primordial scale-invariant power spectrum, the presence of the scalar modes with spectral index Us ~ 1 
and a standard ionization history of the universe. All the power spectra are normalized to COBE/DMR four- Year 
data (Bunn & White 1997). Panels al) and a2) present the CMB anisotropy and polarization power spectra for 
sin^ 2^0 — 0.8, A^j^ = S,^^ — £^1, =4 [L,^ = 3.54) and some values of Am^. Panels bl) and b2) show the same power 
spectra for Am^ = 0.24eV^, A^^ = 4, and some values of sin^ 2f?o- Panels cl) and c2) present the CMB anisotropy 
and polarization power spectra for Am^ = 0.24eV2, sin^ 20o = 0.45 and A^,,/L^ = 0/2.03, 1/2.13, 3/2.88, 4/3.54. 
For all cases we indicate the corresponding values of Vti, (equation 17). Without oscillations our power spectra are 
in agreement with the CMB anisotropy power spectra presented by Lesgourgues & Pastor 1999. Also, in panels bl) 
and b2) we present (dashed lines) the CMB angular power spectra obtained for: fi^ k, 0.015, ^^ = 0, fif, = 0.023, 
r^A ~ 0.7, ilc — 0.262, /iQ — 0.65, one massless and two massive neutrino species. We check that our results match 
exactly those obtained by using the unmodified version of the CMBFAST code. 
The power spectra presented in Figure 1 show two distinct features: a vertical shift of the C; at large I (present in 



L. Popa et al.: On the detection of neutrino oscillations with Planck surveyor 



Fig. 2. The evolution with the scale factor of the energy densities, relative to the critical density, of the various components, 
for the models presented in Figure 1. Panel a): sin^ 25o = 0.8, A^^ — 4 (i^ — 3.54) and Am^: 0.24eV^ (continuous lines), 
0.2eV^ (long-dashed lines), 0.12eV^ (dashed lines), (dott-dashed lines). Panel b): Am^ = 0.24eV^ A^^ = 4 (L^ = 3.54) 
and sin^26lo: 0, 0.23, 0.45, 1 (see also the text). Panel c): Am^ = 0.24, sin^ 2^0 = 0.45 and A^„/L„: 4/3.54 (continuous 
lines), 3/2.88 (long-dashed lines), 1/2.13 (dashed lines), 0/2.03 (dott-dashed lines). Panels al), bl) and cl) show the separate 
contributions from the two types of massive neutrinos (together with antineutrinos) for the same cases of the above panels. 
For all values of the scale factor, ^ f^i ~ 1. 



all panels) with the increasing of Arn^, sin 29o and A^^ values, that results from the modification of neutrino free 
streaming scales, and an horizontal shift of the Doppler peaks to lower I when Arn^ and A^^ are increased, that 
results from the increase of the sound horizon at recombination. The horizontal shift of C; is not present in panels 
bl) and b2) because the variation of sin 26o, with Am^ and A^^ kept fixed, docs not change fl^^ and consequently 
the Hubble expansion rate is not changed, (see equations 14 and 15). 

For a better understanding of these features, Figure 2 shows the evolution with the scale factor of the energy 
densities, relative to the critical density, of the various components {^cdm for CDM, flbar for baryons, ^rad for 
radiation in the form of one massless neutrino and photons, Qi, and Q^^ for the massive neutrinos and antineutrinos, 
flyac for vacuum), for the models presented in Figure 1. Panel a) presents the evolution with the scale factor of 
rii for sin^ 2^0 = 0.8, A^^ — 4 (L,^ = 3.54) and some values of Am"^. The same time evolution is presented in 
panel b) for Arn^ = 0.24eV^, A^^ = 4 {L^ — 3.54) and some values of sin^ 29o. In these cases, as consequence of 
the conservation of the Hubble expansion rate for all values of the scale factor, the time evolution of fli do not 
changes with the variation of sin^ 29q. Panel c) presents the evolution with the scale factor of il^ for Arn^ = 0.24, 
sin^ 29n = 0.45 and some values of A^^/L^. Panels al), bl) and cl) show the separate contributions from the two 
types of massive neutrinos for the same cases of the above panels. For all values of the scale factor, ^^ fli — 1. 

Figure 3 shows the time evolution (in synchronous gauge) of the perturbed energy density of the density field 
components. The left panels of Figure 3 present the evolution with the scale factor of the perturbed energy density 
of different density field components in the presence of neutrino oscillations {Scdm for CDM, dbar for baryons, drad 
for radiation in the form of massless neutrino and photons, di, for v^ plus P^ and S^^ for Vt plus Pt), for the mode 
k « 0.2Mpc^^. The right panels of Figure 3 present, for the same mode k, the evolution with the scale factor of 
the quantity e — cr^j — ct^j (proportional with the variation of the massive neutrinos shear stress) , where ai,-^ and 
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(j^o are: 



^ Pi + ^i . Pj + -Tj 

Here pi, and Pjy are the ncutrino/antineutrino energy density and pressure and Sp^, is the perturbation of neu- 
trino/antineutrino energy density. For all cases the total neutrino mass is m^ = 0.6eV, the total neutrino degener- 
acy parameter is ^jy = 5 and the cosmological model is the ACHDM model presented in Figure 1. In panels al) and 
a2) of Figure 3 sin^ 29a = 0.8 and A^^ ~ 4 are fixed, and Am^ — 0.24eV^ (continuous line) and Am^ = 0.12eV^ 
(dashed line). In panels bl) and b2) of the same figure Am^ = 0.12eV^ and A^^ ~ 4 are fixed and sin^ 29o = 0.8 
(continuous line), sin^ 20o = 0.23 (dashed line). We observe that the growth of 5^. is suppressed when Am^ and 
sin 26q values are increased, while the variation of neutrino shear stress gets larger with the increasing of these 
parameters. In panels cl) and c2) of Figure 3 Am^ — 0, sin^ 29o = and A^^ — 4 (continuous line) and Arn^ ~ 0, 
sin^ 2^0 = and A^j^ — (dashed line). For this last case, Spu ~ 8pv^ and |e| < 10~^. 



Fig. 3. The evolution with the scale factor (in synchronous gauge) of the density perturbations (left panels) and of 
e = CTyj — o"i^2 (sse also the text) for the mode k ~ 0.2Mpc~^. Panels al) and a2): sin^ 26o = 0.8, A^^ — 4 and Am^ = 0.24eV^ 
(continuous line) and Am^ = 0.12eV^ (dashed line); Panels bl) and b2): Am^ = 0.12eV^, A^^ = 4 and sin^ 2610 = 0.8 (con- 
tinuous line) and sin 26^0 ~ 0.23 (dashed line); Panels cl) and c2): Am — 0, sin 26o — 0,A^^ = 4 (continuous line) and 
Am^ = 0, sin'^ 26'o = 0, A^i, = (dashed line). For all cases the total neutrino mass is m^ = 0.6eV, the total neutrino 
degeneracy parameter is ^i, = 5 and the cosmological model is the ACHDM model presented in Figure 1. 



4. Detection of neutrino oscillations with CMB experiments 

Different methods can be used to quantify the performance of a CMB experiment in measuring a given set of 
cosmological parameters. We estimate the errors on the oscillation parameters Am"^ and sin^ 29q and the lepton 
asymmetry L^, by using the Fisher information matrix approximation, widely employed in the literature (e.g. 
Efstathiou & Bond 1999, Popa et al. 1999). 
The Fisher information matrix elements Fij measure the width and the shape of the likelihood function around its 



L. Popa et al.: On the detection of neutrino oscillations with Planck surveyor 9 

maximum (Efstathiou & Bond 1999). The minimum error that can be obtained on a parameter Si when we need 
to determine all parameters jointly, is given by: 



SS^ = VJ^7\ (19) 

depending on the experimental parameter data set and the target model and the class of considered cosmological 

models. 

If only the temperature anisotropy power spectrum Cti is used, the Fisher information matrix reads as (Efstathiou 

& Bond 1999): 

^.-j:'-^-^o.-HCh).'-^. (20) 



I 



Ssi ' '" 5s 



If both anisotropy and polarization power spectra are used, the Fisher information matrix is given by (Zaldarriaga 
& Seljak 1997a; Zaldarriaga 1997): 

F.-j:j:'-§^cov-\c.,cyo'-^, (21) 



I X,Y 



Si OS 



J 



where: X and Y stands for T, E, C and B power spectra and C'ov ^ is the inverse of the covariance matrix. For 
the purpose of this work we assume only scalar modes, then the relevant covariance matrix elements in equations 
(21) ~ (23) are: 

(Zi + ij/sfcy 

C'o^'(C'l;) = ,^, ,^,, (CEi+Wp'Bpf)^, 

Cov{Cli) = -—^-—\Cl, + (Cti + Wp^Bpl){CEi + ^p'Spf )], 



(2/ + \)h^y 

2 


[21 + \)Uky 

2 


[21 + \)Uky 

2 


(2; + \)]sky 

2 


(2; + \)u^v 

2 



-a 



Cov{CtiCei) = ,^, , ^,, C'cu 

Cov{CtiCci) = TT^r^^TTT-^ci^'^Ti+w^^Bpl), 

Cov{CeiCci) = ,^,~, CciiCEi+Wp'Bpf). (22) 

(21 + IJ/sfcy 

Here we set wt — X^c ^c" lo'^ anisotropy and wp — ^^ w^ for polarization (with the sum performed over 
detector channels), w]^ ~ (cj^cpia;)"^ and w^ — (<7c^c,pix)^'^, where crj and a^ are the relative noise per pixel 
for anisotropy and polarization channels (Knox 1995); B^i — ^^ B'^iwJ /wt and Bpi = ^^ B^^w^ /wp account for 
the beam smearing, B^i — e^'('+^)/'s is the Gaussian beam profile. Is — \/81n2(0c)7^;im and fsky is the fraction of 
the sky used in the analysis. 

Assessing Planck performances under realistic assumptions is a very delicate task: both instrumental system- 
atic effects (e.g. Delabrouillc 1998, De Maagt et al. 1998, Maino et al. 1999, Burigana et al. 1998 and references 
therein) and astrophysical contamination (e.g. Toffolatti et al. 1998,1999, De Zotti et al. 1999a, b and references 
therein) have to be carefully understood, reduced by optimizing the telescope (Mandolesi et al. 1998b) and the 
instrumental design (Bersanelli & Mandolesi 1998, Lamarre et al. 1998) and accurately subtracted in the data anal- 
ysis. For simplicity, we consider only the MAP and Planck "cosmological windows" , i.e. those frequencies that 
are minimally affected by foregrounds contaminations, and neglect the degradation in CMB power spectrum recov- 
ering introduced by foreground contaminations and instrumental effects. As shown in Figure 4, in the cosmological 
channels extragalactic source fluctuations are comparable to CMB ones only at I > 10^, both for temperature (left 
panel) and polarization (right panel) anisotropics, whereas galactic fluctuations, which power spectrum overwhelms 
the extragalactic foreground one at / ^ 10^, are below the CMB ones at least far from the galactic plane. 

Then, the possible imprinting of neutrino oscillations in the CMB anisotropy are not significantly masked from 
astrophysical contaminations at frequencies ~ 70 -^ 200 GHz. We exploit then separately different sets of Planck 
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data: the "cosmological" LFI channels alone (70 and 100 GHz channels), the "cosmological" HFI channels alone (143 
and 217 GHz channels and 100 GHz but for the temperature fluctuations only) and both together. For comparison, 
we consider also the set of data from the two MAP "cosmological" channels at 60 and 90 GHz. Table 1 lists the 
experimental parameters of the various experiments that we considered in our calculations. 



Fig. 4. Temperature (left panel) and polarization (right panel) power spectra for the neutrino oscillation target model 
considered here (see also the text) compared to astrophysical component power spectra: the galactic temperature fluctuations 
at moderate galactic latitudes, sum of synchrotron, free-free and dust fluctuations, and the extragalactic source temperature 
fluctuations modelled as in Toffolatti et al. (1998) and De Zotti et al. (1999a) and the polarization fluctuations from galactic 
synchrotron and dust emission and separately from extragalactic radiosources and infrared sources modelled according to 
De Zotti et al. (1999b) and references therein. We report also the white noise nominal power spectrum for LFI and MAP, 
taking into account the effect of resolution loss at high I due to the beam convolution. 

The Planck and MAP data will be mainly used for determining the cosmological parameters; we exploit then 
the Fisher matrix method by evaluating together the uncertainties {la errors) of a simple set of cosmological 
parameters and of the parameters characterizing the neutrino oscillation model discussed here. As well known, the 
errors quoted for the considered parameters depend on the experiment sensitivity but also in part on the assumed 
target model and the set of considered parameters. 

We assume as target model a spatially flat ACHDM model without oscillations having: flf, = 0.023, Ha ~ 0.7, 
D,c ^ I - rib - ^u - ^A, ho = 0.65, Us ^ I, rriy ^ m^^ + m^^ = 0.6eV, ^^ = ^^^ + ^^^ = 5, A^^, = 1 (L^ = 2.13), 
Arn^ ~ 1.24 X lO^^eV^, and sin^ 29o = 0. We also assume a scale invariant power spectrum, the presence only of the 
scalar modes having a spectral index Us — 1, and ignore the contribution of the tensorial modes and reionisation 
effects. The tensorial modes have a small contribution only to the polarization signal at low order multipoles 
(Crittenden & Turok 1995, Seljak 1997, Kamionkowski & Kosowsky 1998) that will be difflcult to detect against 
the polarized foregrounds. Also, the reionisation effects can help in breaking the degeneracy only in the presence of 
the polarization signal for a large value of the optical depth to Thomson scattering (Zaldarriaga, Spergel & Seljak 
1997). 

We numerically compute the partial derivatives of the power spectra with respect to each relevant parameter Si 
by approximating them with the (possibly central) finite differences between the two power spectra obtained by 
varying with a quantity zLAsi/2 the parameter Si in consideration around its value Si^ assumed in the considered 
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target model. Typically, Asi is taken in the range between few percent, in order to have a step small enough for an 
accurate evaluation of the derivative and large enough for having quite sensitive and numerically stable variations 
of the power spectrum. 

Table 2 presents 1-a errors on the estimates of the relevant cosmological parameters and neutrino oscillation 
parameters obtained from the anisotropy alone and anisotropy plus polarization, for the experimental parameters 
listed in Table 1 and few values of fsky ■ The errors presented in Table 2 arise because of the geometrical degeneracy 
(see Efstathiou & Bond 1999 and the references therein) that imposes limits on the determination of ho and J^a, 
and the correlations existing among the energy density parameters and Us and sin 26o. 

These errors can be understood in terms of the ability of each experiment to determine the properties of Doppler 
peaks (Efstathiou & Bond 1999): the location of the first Doppler peak and the positions, heights and relative 
amplitudes of the first and subsidiary Doppler peaks. The Doppler peak positions Id [Id ~ fn'^dA{z)/rs, where 
dA{z) is the angular diameter distance to the last scattering, r^ is the sound horizon distance and m is the Doppler 
peak number] suffer the degeneracy between Qa and ho when Qm is fixed {Qm = ^t + ^c + ^i^ = 0.3). From 
equation (19) 51,y value is given by the Atti^ and Li, values when the total neutrino mass is fixed, and therefore 
Ato^ and L^, are also correlated. The Doppler peak heights suffer the degeneracy between flc, ^b, ^ly, ns and 
sin^26'o. 

Clearly, the Planck LFI, probing the subsidiary Doppler peaks structure to high multipoles / ^ 1500, will produce 
a significant improvement with respect to MAP (sensitive up to Z ^ 1000) in the determination of all parameters 
reducing their uncertainties. A further improvement can be achieved by exploiting the better sensitivity and angular 
resolution (up to / ^ 2000) of Planck HFI, the full Planck performance at very high / being limitated essentially 
by the foreground contamination. 

On the other hand, other kinds of astronomical observations provide information on the cosmological parame- 
ters, complementary to those derived from CMB anisotropies. Galaxy cluster observations limit the shape parameter 
r and the normalization of the galaxy power spectrum ag and, combined with the galaxy peculiar velocity fields, 
provide informations on the matter density in the universe flm (Dekel 1994, Strauss & Willick 1995). Weak gravita- 
tional lensing provides also direct estimates of qo (see e.g Seljak 1998). The ratio between baryon and total mass in 
clusters flh/flm ~ (0.01 - 0.02) ho"^ (see e.g. Dodelson et al. 1996 and the references therein), informations from the 
BBN [Qb = 0.02/i^^ (Buries & Tytler 1998)] and the Hehum Lyman- Alpha Forest (Wadsley et al. 1999) constrain 
flm or ho if ^m is known from other kinds of measurements. The complementarity between Type la supernovae 
luminosity distances and CMB anisotropy provides a powerful tool in breaking the degeneracy between flk and JIa 
(see e.g. Perlmutter et al. 1998, Efstathiou & Bond 1999, Efstathiou 1999 and the references therein). The ages 
of the oldest globular clusters set lower limits on the age of the universe (see e.g. Chaboyer et al. 1998) and then 
on Qk and J^a if ^m is well constrained. For a spatially flat universe this contributes to break the geometrical 
degeneracy in ho - ^a plane. For CDM models with known baryonic content, measurements of F from large scale 
structure, combined with independent determinations of firm provide estimates on ho (Efstathiou & Bond 1999). 
Constraints on the scalar spectral index n^ can be inferred also from limits on the primordial black hole abundance 
(Green, Liddle, & Riotto 1997) and on the CMB spectral distortions (Hu, Scott, & Silk 1994). 

Recent CMB anisotropy measurements from Boomerang and MAXIMA-1 have been in fact used jointly to 
other astronomical observations to determine cosmological parameters through the so-called analysis with prior 
(Lange et al. 2000, Balbi et al. 2000). Of course, adding informations from other kind of observations helps both 
in circumventing the degeneracy problem and in improving the accuracy in the determination of the cosmological 
parameters. 

Such a kind of analysis is out of the aim of this work. We just estimate here lower limits on the MAP and 
Planck sensitivity in determining the neutrino oscillation parameters, by exploiting the Fisher matrix method, 
assuming that all the cosmological parameters are known and estimating only the joint errors on Arn^, sin 26o, 
and Li, that can be obtained from CMB anisotropy plus polarization data. 

Table 3 presents these errors {la errors) for Planck (LFI and HFI together) and MAP experiments; together 
with the results shown in Table 2, it provides a range of Planck and MAP sensitivity to neutrino oscillation 
parameters. Of course, in this case we find a significant reduction of the errors on the neutrino oscillation parameters 
for both Planck and MAP. These values can be also considered as estimates of lower limits on the errors on the 
neutrino oscillation parameters, as they could be in principle stetted by Planck and MAP experiments. 

Figure 5 presents few confidence regions of the neutrino oscillation parameters space that can be potentially 
detected by the Planck surveyor and few allowed regions of Am^ and sin^ 29o at the same confidence level (CL) 
obtained by MACRO CoUaboration (Ambrosio et al. 1998). 
According to the standard x^ method (sec e.g Particle Data Group 1998) an input model defined by the parameters 
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(Ato^, sin^ 29o) is more likely to have occurred higher is the probability to obtain x^ values larger then a specific 
value Xo- The cumulative probability that defines a certain confidence region on the parameters is given by: 

CLixl)^ i(x^A™^sin2 20o)dx^ (23) 

Jxl 

where Xo — Xmin + ^X^) "^ith ^X^ — 4.61, 9.21 at 90% and 99% CL respectively (for a x^ distribution with two 
degrees of freedom). The likelihood function in the equation (25) is defined as (Efstathiou & Bond 1999): 

L \ y {Ci(s) - Q{sO))^ 

i<i„ 



.21n^^ = y 



where C((so) is the power spectrum for the target model, Ci{s) is the power spectrum for the input model having 

the same cosmological parameters as the target model and Am^ and sin^ 20o are in the intervals (0 — 0.36)eV^ and 

(0—1) respectively. 

Panel a) of Figure 5 presents the confidence regions that could be obtained with the Planck surveyor for different 

fractions of sky. The target model is the KCHDM model used in the previous section. Panel b) of the same figure 

presents the confidence regions obtained with Planck for fsky — 0.65 and the same target model and for another 

target model. 

In each panel we also present the allowed regions of Atti^ and sin^ 26*0 as derived from the MACRO experiment. 

We conclude that exists a significant overlap between the region of the oscillation parameter space that can be 

potentially detected by Planck surveyor and that implied by the atmospheric neutrino oscillations data. 



Fig. 5. Confidence regions of the neutrino oscillation parameter space that can be potentially detected by Planck surveyor 
by using CMB anisotropy measurements. Panel a): the allowed confidence regions at 99% CL (A) and at 90% CL (B) 
obtained for f^ky = 0.65 (continuous lines) and fsky ~ 0.8 (dashed lines). The target model has Am,^ = 1.24 x 10~'^eV'^ and 
sin 26o = (see also the text). Panel b): the allowed confidence regions at 99% CL (A) and at 90% CL (B) for fsky ~ 0.65. 
The target model has: Am^ = 1.24 x 10"^eV^ and sin^ 20o = (continuous lines) and Am^ ^ lO'^^eV^ and sin^ 26*0 = 
(dashed lines). In each panel we also present (dash-dotted lines) the allowed regions of Am^ and sin^ 2So at 99% CL (Al) 
and 90% CL (Bl) measured by the MACRO experiment. 
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5. Conclusions 

The imprint of relic degenerated neutrino oscillations on the Cosmic Microwave Background (CMB) angular power 

spectra is evaluated in a ACHDM model consistent with the structure formation theories, allowing in the same 

time a pattern of neutrino masses consistent with the atmospheric neutrino oscillations data. 

Under the assumption that the total neutrino mass is known, as derived by the galaxy redshift surveys (Hu et 

al. 1998), we show that oscillations occurring between two active degenerated neutrino flavors leave detectable 

imprints on the CMB anisotropy and polarization power spectra. 

By using the CMB anisotropy measurements, we find that the relic degenerated neutrino oscillations can be detected 

by Planck surveyor at la level if the present value of the lepton asymmetry is L^ > 9.7 x 10~^, the difference 

of the neutrino squared masses is Am^ > 2.91 x lO^^eV'^ and the vacuum mixing angle is sin^ 29o > 0.46, when 

we consider other cosmological parameters that simultaneously affect the CMB angular power spectra and a sky 

coverage of fsky = 0.8. A slight improvement can be obtained by using the CMB anisotropy and polarization data. 

Assuming that all the other cosmological parameters may be known from other kind of astronomical observa- 
tions, we estimate lower limits on future CMB anisotropy and polarization experiment sensitivity in determining 
the neutrino oscillation parameters alone: as shown in Table 3, the errors on neutrino oscillation parameters de- 
crease at least by a factor two with respect to the case in which no prior on cosmological parameters are assumed. 
These values provide estimates of lower limits on the errors on the neutrino oscillation parameters, as they could 
be in principle stetted by Planck and MAP experiments. 

Finally, we compare the confidence regions of the neutrino oscillation parameter space that can be potentially 
detected by the Planck surveyor with the few allowed regions of Arn^ and sin^ 26*0 obtained by MACRO Collab- 
oration (Ambrosio et al. 1998). We find a significant overlap between the region of the oscillation parameter space 
that can be potentially detected by Planck . 
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Table 1. Summary of instrumental performances of future CMB space missions in their "cosmological channels". For MAP 
and Planck LFI we compute af = \/2o-J (Zaldarriaga & Soljak 1997). 



Instrument 


u{GRz) 


dfwhm 


a^/i),K/K) 


(™^)-Vio-^^ 


a^/{fiK/K) 


(»r)-vio-^^ 


MAP 
(Bennet et al. 1996b) 


60 
90 


21' 
12.6' 


12.1 
25.5 


5.4 
8.7 


17.11 
36.06 


10.92 
17.46 


Planck LFI 
(Mandolesi et al. 1998) 


70 
100 


14' 
10' 


3.6 
4.3 


0.215 
0.156 


5.09 
6.08 


0.429 
0.312 


Planck HFI 
(Puget et al. 1998) 


100 
143 
217 


10.7' 

8' 

5.5' 


1.7 
2.0 
4.3 


0.028 
0.022 
0.047 


3.7 

8.9 


0.074 
0.020 



Table 2. la errors on the estimates of the cosmological parameters and neutrino oscillation parameters obtained from the 
power spectrum statistics. 
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MAP 


Sho/ho 


0.11 


0.13 


0.14 


0.10 


0.12 
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1.95 


2.16 


0.64 


0.72 


0.79 




SL^/L^ X 10^ 


4.09 


4.57 


5.07 


1.63 


1.84 


2.04 




SAm^/Am^ 


2.11 


2.35 


2.61 


1.39 


1.56 


1.73 




(5sin2 26io 


0.41 


0.46 


0.51 


0.37 


0.41 


0.46 



Table 3. la errors on the estimates of the neutrino oscillation parameters and lepton asymmetry obtained from anisotropy 
and polarization power spectrum statistics. 



J sky 




1. 


0.8 


0.65 




5L^/L^ X 10'' 


1.21 


1.35 


1.50 


MAP 


SAra^ /An? 


4.08 


4.56 


5.06 




S%m^2eo 


0.68 


0.77 


0.85 




SL^/L^ X 10* 


2.13 


2.38 


2.64 


Planck 


SAn? /An? 


0.79 


0.88 


0.97 




8 sin^ Wo 


0.13 


0.19 


0.21 



